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Abstract 

We study the evolution of an oscillator interacting via the most general bilinear coupling (with 
time-independent coefficients) with an "environment" consisting of a set of other harmonic oscilla- 
tors. We are mainly interested in a possibility of using the Fokker-Planck equation to describe this 
evolution. Studying different interaction Hamiltonians, we show that unambiguous reduction to the 
Fokker-Planck equation is possible only within the framework of the so called rotating-wave approxi- 
mation. As special cases we consider in detail the evolution of two coupled oscillators and relaxation 
of a charged oscillator in a uniform magnetic field. 
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1 Statement of the problem 

This part is devoted to the multidimensional generalizations of the harmonic oscillator model. As a rule, 
these are various systems of interacting oscillators, some of which may be placed in external uniform 
electric and magnetic fields. More precisely this class of systems can be characterized by a Hamiltonian 
that is assumed to be in the form of the general inhomogencous quadratic form 

<N 

£=-qBq+Cq. (1.1) 

& : 2 

Here, q is conceived as a 27V-dimensional vector (TV being the number of degrees of freedom) whose 
components are linear combinations of Cartesian coordinates and momenta conjugated to them. The most 
frequent choice is q = (pi, . . . ,pn, %i, ■ ■ • , a; at), but in the presence of an external magnetic field it is more 
convenient to deal not with the canonical momenta pj but with the kinetic momenta Hj = pj — eAj (x) /c, 
where A(x) is the vector potential. Other choices are also possible, taking into account the concrete 
physical applications. For instance, the components of the vector q may be bosonic annihilation and 
creation operators constructed from the coordinate and momentum operators, etc. B is a symmetric 
27V x 27V matrix, which may depend on time, as also the 27V-dimensional vector C. 



X 

Multidimensional systems with Hamiltonian (1.1) were the subject of investigation in numerous pa- 
pers: see, e.g., |?], ||, Q-|l^| and the references therein. In the present paper we consider the case where 
the multidimensional vector q may be split in two parts: q = (Q,£), where the vector Q describes a 
small subsystem, whereas the vector £ is related to a "thermostat" . In the simplest cases the evolution 
of a unidimensional harmonic oscillator coupled with a "thermostat" was studied, e.g., in |]|,(l7|]-|3§]. 
We investigate the most general quadratic "interaction Hamiltonians" and compare the results of this 
"microscopic" approach with different phenomenological models considered in reviews Emphasis 
is placed on application to the problem of harmonic oscillator relaxation, including the case where a 
(charged) oscillator is placed in a uniform magnetic field. 

2 Phenomenological Fokker-Planck equation 

The evolution of any closed system is governed by the quantum Liouville equation for the statistical 
operator p, 

ihdp/dt = Hp-pH. (2.1) 



For any quadratic Hamiltonian (1.1) this equation results in the linear equation for the average values of 
the 27V-vector q: 

(q) = -SB(q) - EC, (2.2) 
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where (q) = Tr(j5q), and antisymmetric nondegenerate 2N x 2N matrix £ with c-number coefficients is 
defined via the commutation relations between the operators q a , 



a,/?= 1,2, 



,2N. 



(2.3) 



The operator equation (2.1) can be transformed into a partial differential equation upon the choice of 
some concrete representation for the statistical operator. For instance, in the coordinate representation 
we write q = (p,x), p and x being the N- vectors, and split the matrix B into N x N blocks: 



bi b 2 
b 3 b 4 



bi =b x 



b 2 = b 3 



b ■ — \\b mn \ 



(the tilde designates a transposed matrix). Then we obtain the following second-order equation for the 
density matrix p(x,x'): 



. % dp 



d 2 P 
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d 2 p 



ihb n 2 in 



dp , dp 



' dx r . 



■dx> 



+ (XrnXn ~ %' m X 'n) P ~ ^/OTrb 2 



\ ^ ^ rn 



dp_ 

dx' 

w ^ rn 



+ c 2 ( x m x rn) ' 



(2.4) 



Here Ci and c 2 are A^-dimcnsional components of the vector C = (ci, c 2 ). The disadvantage of Eq. (2.4) 
is the broken symmetry between the coordinates and momenta, which is inherent in the Hamiltonian. 
This symmetry can be restored if one proceeds from the complex density matrix to the real Wigner 
function 



W(p,x)= j p(x + £/2,x-£/2)exp(-ip£//i)d& 
p(x,x') = j w(p^(x + x'))exp[ip( X - X ')/h} dp/(2nh) N . 



(2.5) 
(2.6) 



Applying transformation (2.5) to Eq. (2.4) we obtain an equivalent equation that is first order with 
respect to all the derivatives: 



dW 



dW 

dp m 



(Pnb? 



+ cT) 



dW 

dXm. 



(2.7) 



This equation demonstrates the distinction and the advantage of the Wigner function for the descrip- 
tion of quadratic quantum systems (other remarkable features of the Wigner function were discussed, 
e.g., in [Eol). Eq. (2.7) assumes an especially compact form in terms of the 2A^-vector q, 



dW 

~~dT 



d 
dq a 



[(EBq + EC) Q W] 



Here the 27V x 2N matrix S, in accordance with Eq. (2.3), equals 

E 






In 


-In 






(2.8) 



(2.9) 



In being the N x N unit matrix. 

Now let us suppose that we have a given linear equation for the first-order average values 

(q)=A(t)(q)+K(t), (2.10) 

with an arbitrary matrix A(t) and an arbitrary vector K(t). The problem investigated in this section is 
whether it is possible to find an equation for the statistical o perat or or the Wigner function that would 
result in the given Eq. ( 2.10 ). It is trivial to check that Eq. ( |2.10|) is the consequence of the equation 



dW d ,, . _ _v Trrl 

lH=-W a [iM + K)aW] 



(2.11) 
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However, although Eq. ( [2.11 ) preserves the normalization of the Wigner function, 



W(q) dq/(2irh) 



N 



(2.12) 



it does not agree, in general, with the fundamental quantum mechanical principle of the positive definite- 
ness of the statistical operator and (this is almost the same) with the uncertainty relations. Indeed, let 
us consider the operator F = otj (qj — (qj)) with arbitrary complex coefficients ctj. For any stattistical 
operator, due to its nonnegative definiteness, the inequality (F'F) = Tr (pF'F) > must hold. Taking 



into account the structure of the operator F, we arrive at the conclusion on the nonnegative definiteness 
of the bilinear Hermitian form 



®jk = ((% ~ (qj)) (& - <»»)• 



whose matrix $ is constructed from the centered second-order moments of the operators qj , j,k 
1.2.... , 2N. It is convenient to distinguish even and odd parts of the matrix $. The symmetric part 
consists of the symmetrized second moments (covariances) 



Mij = Mji = ~ (%q k + %%) - (<lj){<lk)- 



(2.13) 



The antisymmetric part, in accordance with Eq. (2.3), is expressed through the commutator matrix, so 
that 

The elements of the matrix M. are calculated in terms of the Wigner function as follows: 



Mi 



q i q j W(c l )dc l /(2nh) 



N 



qW(q) dq/(27r7i) 



iV 



Due to Eq. ( 2.11 ), the variance matrix satisfies the equation 



(2.14) 



is positive (the simplest 



M = AM+MA. 

The first-order averages do not influence the variances in the case under study. 
If the coefficients ctj are chosen in such a way that the commutator F, F^ 

example is F = x + ip) , then the operator F coincides within a constant factor with a boson annihilation 
operator. Choosing the initial state to be the vacuum state for this operator, at the initial instant we 
have 

Tr (F^Fp(ofj = a*$(0)a = 0. 
Moreover, the matrix equalities <!>(0)q; = a*$(0) = hold as well. For the chosen initial state we have at 



t > 0, due to Eq. (2.14) 



a*<f>(t)a 



a 
1 



M(t)-yE 



A(0)7W(0)+7W(0)A(0) ta + 0(t 



-ihta* 



A(0)E + SA(0) a + 0{t 2 ) 



(2.15) 



For the Hamiltonian systems we have A = — EB, A = BE, and the linear with respect to time term 
disappears. But for an arbitrary matrix A the right-hand side of Eq. ( p. 15 ) can be negative. For example, 
in the model of a damped oscillator with the equations for the averages 



x = p, p= -uj x - 2jp, 



(2.16) 
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the matrices entering Eq. ( 2.15 ) read 



-27 







AE + EA 




2- 



-27 




(2.17) 



Then, for the initial vacuum state of operator F = x + ip, when a = («, 1), we obtain from Eq. (2.15) 

a*${t)a = -2^7 + 0(t 2 ) < 0. 



Consequently, in the general case Eq. ( p. 11 ) is unacceptable. Thus we need more complicated gener- 
alizations. The simplest possibility is to add terms with second derivatives to the right-hand side of Eq. 



(2.11), i.e., to transform this equation into the Fokker-Planck equation: 



at Oq a oq a dq f 3 



(2.18) 



where the diffusion coefficients D a p — Dp ai combined into a symmetric matrix D = ||_D a ^||, may depend 
on time but do not depend on the coordinates. The new "diffusion" terms do not change the equation 
for the average values ( |2.10 ), moreover, they preserve the normalization ( 2.12 ). But they enable one to 
"save" the nonnegative defmiteness of the statisti cal o perator. Indeed, considering the evolution of the 
bilinear form a*$(t)a, we obtain, instead of Eq. (2.15), the equation 



a*<S>(t)a = ta* [2D + — A(0)E + £A(0) 



0(t 2 ), 



since Eq. ( [2.14 ) is replaced by 



M = AM + MA + 2D. 



(2.19) 



Consequently, the necessary condition of the compatibility of Eq. (2.15) with the principles of quantum 
mechanics is the nonnegative defmiteness of the matrix 



D* = D 



ih 
T 



AE + EA 



> 



(2.20) 



at any instant of time. Moreover, it can be proved |39| |41J that the condition D* > is sufficient as well. 

Vector q was defined above as q = (p, x), and the matrix £ had the explicit form given by Eq. Q2.9|) , 
Let us make the time-independent transformation of the variables 

q' = Tq, det T^O, ImT ^ 0. 



Then Eqs. (fc.lOD and (|2.18|) preserve their forms, provided the matrices A and D are replaced by the 
matrices A' = TAT -1 and D' = TDT. The nonnegative defmiteness of the matrix D* of Eq. (2.20) 
is equivalent to the nonnegative defmiteness of the matrix D'„ = TD*T. The latter has, in turn, again 



the form ( 2.20 ), if one replaces the matr ix E by £' = TET, but this is just the transformation law of 
any matrix defined according to Eq. (2.3). This way we arrive at the important conclusion, that all the 
formulas containing the matrices A, D, and E are valid not only in the case where the components of 
the vector q coincide with the canonically conjugate momenta and Cartesian coordinates, but also in 
the general case where the components of the vector q are arbitrary Hermitian operators with c-number 
commutators, provided the matrix E is defined according to Eq. (|2.3| ). 

To transform Eq. (2.18) to an operator form that is independent on the concrete representation, one 
should take into account the following correspondence relations between the operators qp, pq, and their 



Weyl symbols (they result from Eqs. (2.5), (2.£)): 



qp ■ 



pq 



qW(q) i (qp + pq), ^ ^E 1 (qp - pq) 



(2.21) 
(2.22) 



Making transformations ( 2.22 ) in Eq. ( ^.18| ) we arrive at the equation 
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dp 

dt 



I 

i 



qE 1 Aqp + pqAE _1 q + q ^E _1 A + AS" 1 ) pq 
[qE^Kp - jSqE^K] - -j [qSqp + ,SqSq - 2qSpq] 



where the symmetrical matrix 



S = -E _1 DE 



must satisfy a constraint equivalent to Eq. (2.20): 



S* = S — — 



AE~ 



E _1 A 



> 0. 



(2.23) 
(2.24) 

(2.25) 



Eqs. ( 2.21 )-( 2.25 ) hold for any vector q w ith a c- number commutator matrix E of (2.S). 

Comparison of the elegant equation (2.18) for the Wigner function with the much more cumber- 
some equation (2.23) demonstrates once more the advantage of the Wigner representation for describing 
quantum systems with linear equations of motion for the averages. 

There is an important difference between the quantum Fokker-Planck equation (2.18) and its classical 
counterpart. The classical Fokker-Planck equation contains, as a rule, second derivatives only with respect 
to momenta. However, such a simple set of the diffusion coefficients is unacceptable in the quantum case 
(although sometimes this incorrect equation was considered: see, e.g., p2[ ). This statement can be easily 
demonstrated on the example of system (2.16). Writing the matrix D in the form 



D 



D„ D„ 



D, 



and taking into account Eq. (2.17) we obtain the matrix 



D, 



D 



p 

ih-y/2 



D 



P x - ihj/2 



The condition of its positive definiteness is given by the inequality 



det D» = D P D X - Dp X - h 2 j 2 /4 > 0, 



whose violation lead s to the violation of the uncertainty relations [j43| , |47|| . 
dimensional system ( 2.10 ) the condition D* > is equivalent to the inequality (4 



(2.26) 

For an arbitrary one- 



detD > ?i 2 (TrA) 2 /i6. 



(2.27) 



3 Fokker-Planck equation for a subsystem 

In the preceding section we have shown that any given equation ( |2.10D may be considered as a consequence 
of some suitable Fokker-Planck equation for the Wigner function. The only problem is to select the set of 
diffusion coefficients satisfying the condition D» > 0. It is clear that by taking sufficiently large diffusion 
coefficients one can always satisfy this condition. The problem of finding "minimal admissible" diffusion 
coefficients for some simple systems (such as a one-dimensional harmonic oscillator or a two-dimensional 
isotropic oscillator in a uniform magnetic field) was investigated in Q . Q , Q . 

Here we investigate the following problem. Suppose we have a large closed quantum system with TV 



degrees of freedom, described by Hamiltonian (1.1). Let us split the vector q in two parts: q = (Q,£), 
where the n-dimensional vector Q describes a subsystem, while the vector £ relates to a reservoir. The 
question is: what kind of equation describes the evolution of the subsystem if one performs an averaging 
over the variables of the reservoir? 

The Wigner function of the whole system is given by the relation 

W(q,t)= J G(q,q',t)W(q',0)dq', (3.1) 
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where the propagator G(q, q', f) satisfies Eq. (2.8) and the initial condition G(q, q', 0) = <5(q— q'). Since 
Eq. ( |2.8| ) is first-order with respect to all the variables, its propagator is extremely simple: 

G(q,q',t) = *(q-q.(*;q')), (3-2) 

where the vector 

q,(t;q')=R(i)[q'-A(t)] (3.3) 

is the solution to the classical equation of motion ( |2.2[ ), satisfying the initial condition q*(0;q') = q'. 
Consequently, the 2N x 2N matrix R(i) satisfies the equation 

R= EBRee.4R, (3.4) 

and the initial condition R(0) = l2jv- The vector A(t) equals zero at t = 0. For t > it is determined 
from the equations 

A = R X EC = SRC. 



These two forms are equivalent due to the identities 

R(t)E~ 1 R(t) ee E -1 , R(t)ER(t) ee E, 



R 1 ee ERE 



(3.5) 
(3.6) 



which follow from Eq. (3.4). 

It is clear that the Heisenberg equation of motion for the operator q coincides with Eq. ( |2.2| ), 
provided (q) is replaced by q. The solution to this equation is given by Eq. (3.3) with carets over q 
and q'. Therefore, the identities (3.6) mean nothing but the conservation of the commutation relations 
(2.3) in time, or, in other words, the canonicity of transformation ( |3.3| ) and the unitarity of the evolution 
operator. 

Now let us proceed to the averaging of the total Wigner function over the reservoir variables £. Our 
first assumption is that the initial total Wigner function is factorized: 

W(q,0) = Wb(q)Wi(0- (3.7) 
The second assumption concerns the initial Wigner function of the reservoir. We assume it to be Gaussian, 

Wi(0 = n M (detF)- 1 / 2 exp 



-^-^f-^-t) 



(3.8) 



with some symmetric positive definite 2M x 2M matrix F and a 2M-vector 7 (M being the number 
of degrees of freedom of the reservoir). In particular, it may correspond to a mixed equilibrium state 
[|, 16, 5[J or to a pure squeezed coherent state jj], [5l|| . For any physically admissible Gaussian Wigner 
function the covariance matrix F must satisfy a set of conditions expressing generalized uncertainty 
relations. The simplest among them is the inequality [52| detF > (h 2 /A) M . Moreover, the parameter 

H= (?i/2) M (detF)- 1 / 2 < 1 (3.9) 

characterizes "the degree of quantum mechanical purity" of the Gaussian state: 



(3.10) 



/j = Tip 2 = J W 2 (£) d£/(2Trh) A 
To calculate the averaged Wigner function 

Wt(Q,t) = J W(Q,0d£/(2irh) 



(3.11) 



we split the matrix R and the vector A into rectangular blocks in accordance with the decomposition 
q=(Q,£): 



R 



Rn R12 
R21 R22 



A 



A. 



Q 



A, 



(3.12) 



Then Eqs. (|3_l|), (|3_2|), (|3_7|) lead to the integral 

W e (Q,t) = yj(Q-Rn[Q'-A Q ]-R 12 [e'-A € ]) 



xS (f - R21 [Q' - A Q ] - R 22 [e' - A c ]) 
xW (Q / )W 1 (e)dQ'ded(/(2TTh) M 



(3.13) 
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The integration over d£ is trivial: it simply removes the second delta function. To perform the integration 
over d£' we replace the first delta function with its integral representation (recall that the dimension of 
the vector Q is 2n), 



5(x) 



d^k/(27T 



,2/) 



(3.14) 



Thus we obtain two Gaussian integrals (the first one over d£' and the second over dk) that can be 
calculated exactly due to the well-known formula 



J dxexp(-xAx + bx) = [det(A/7r)]- 1/2 exp QbA _1 b^ 



(3.15) 



Finally we arrive at the same equation (3.1), but with the variables Q, Q' instead of q, q', W^(Q,t) 
instead of W(q, t), and with the averaged propagator 



G«(Q,Q',t) 



{2ir)- n [detM4t)r 1/2 
1 



x exp 



(Q - RnQ' - 5*) AC 1 (Q - RnQ' - 8, 



(3.16) 



where the symmetric matrix M*(t) equals 

M*(t) =Ri2(t)PRi2(t). (3.17) 

Recall that the dimension of the rectangular matrix R12 is 2n x 2M. The vector 6*(i) equals 

$*(t) = R12W7 " Rn(t)A Q (<) - R 12 (i)A 5 (i). (3.18) 

Direct inspection shows that the propagator J3.16] ) (and, consequently, its convolution with any initial 
function) satisfies the Fokker-Planck equation ( 2.18 ) with the following drift matrix A and vector K, 



A — RnR 11 1 , 
K = L- AS*. 



The matrix of diffusion coefficients reads 



1 



D = | \M* - AM, - M«A) = sym [(R 12 - R 11 R^ 1 1 R 12 ) FR 12 



where we have introduced the notation 



symA = - ^A + A^ 



(3.19) 
(3.20) 

(3.21) 
(3.22) 



Thus we have proved that the evolution of the Wigner function of any subsystem of a closed quadratic 
quantum system is governed by some effective Fokker-Planck equation, provided the initial state of the 
"remaining part of the system" is Gaussian. However, the coefficients of the equation obtained in this 
way depend, as a rule, on time even for a time-independent Hamiltonian of the closed system. 



4 Two coupled oscillators 

Let us illustrate the formulas of the preceding section, applying them to a very simple special model, 
where both the "subsystem under study" and the "reservoir" are the harmonic oscillators with a single 
degree of freedom. This example admits exact solutions, in contrast to the more realistic situation of a 
reservoir with a very large number of degrees of freedom, where one has to make various simplifications 
and approximations to obtain a closed result that would be easy to analyze. Consider first the free 
Hamiltonian 

TT 2 2 1 1 1 22 (A 1 \ 

Ha= 2^ + 2 mi ^ + 2^ + 2 m « (41) 
and the most general quadratic interaction Hamiltonian 

-^int = 9 PP PiP2 + g P xP\X2 + gx P xiP2 + 9xxXix 2 . (4.2) 
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We assume all the coefficients to be time-independent. In the notation of the preceding section we should 
write 

q = (P1,X1,P2,X 2 ) , Q=(pi,Xl), £, = (P2,X2)- 

Then the matrix A = -SB (see Eq. (O)) reads 



A 



Its characteristic equation turns out to be biquadratic: 








9xp 


9 XX 







9pp 


9px 


9px 


9 XX 







9pp 


9xp 


m 2 1 






det(.4 - iwl) 



/ 9 2 
Wi + Wo 



2A) u z + uj\ujI + t±. A - g = 0, 



where 



det A\ 2 — QppQxx QpxQxpi 
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9 XX 



The solutions of Eq. (| 



m\m,2 
are as follows 



(4.3) 



(4.4) 

(4.5) 
(4.6) 



1 



(^ + c 2 2 ) + (^2 + A 2 -.g) 1/2 + A 



± 



- (w 2 + c^ 2 ) - (w 2 ^ 2 + A 2 - g) 1/2 + A 



(two other solutions are equal to — u>±). 

Thus, the evolution of two coupled harmonic oscillators can be described explicitly for quite arbitrary 
quadratic interaction Hamiltonians with timc-indcpcndcnt coefficients. If both frequences given by Eq. 
(4.7) arc real, then the particles perform harmonic oscillations. For certain parameters complex normal 
frequencies are possible. Then the motion becomes aperiodic. However, since any normal frequency w + or 
cij_ is accompanied by the frequency with the opposite sign, it is impossible to obtain damped oscillations 
of cither particle. The coordinate and momentum of any oscillator will increase with time. 

To illustrate this statement, let us consider first the case (which see ms the most natural) of the 
interaction via the coordinates, where the only nonzero coefficient in Eq. (4.2) is g xx . Then A = 0, and 
for sufficiently strong coupling, when g > <x> 2 w 2 , we have the real frequency w+ and the pure imaginary 
frequency lj- : 

„1 1/2 



1 



iX, 




Solving Eq. (3.4), we obtain the following formula for the matrix Ru(t): 

p+ cos tot + p- cosh Xt mi(—ujp + sin tot + Xp- sinh Xt) 

^-j- (^j- sin tot + £=■ sinh Xt) p + cos ujt + p_ cosh Xt 



Rn 



where 



1 I,, 
P±= 2±4^ 



(4.7) 
(4.8) 

(4.9) 
(4.10) 



We see that the oscillations of the first particle actually increase. Its energy is derived from the nonpositive 
definite potential energy of the whole system, since the whole system turns out to be unstable when 
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g > ijj\uj\. Note that E g. (|4.9| ) still holds even if g xp ^ 0. If g xx — g xp = 0, but g pp ^ and g px ^ 0, we 
obtain, instead of Eq. (4.9), the expression 



Rn = 



p + cos tot + p- cosh Xt 
j^jj (ujp + sin tot + Xp- sinh Xt) 



-m\Lo\ (^j- sin tot + ^- sinh At) 
p+ cos + p- cosh Xt 



(4.H) 



with the same values of u> and A (provided, of course, g > io\io?£). 

Both frequencies uj + and io_ have nonzero imaginary parts provided the argument of the first square 
bracket in Eq. (4.7) is positive while the argument of the second square bracket is negative. Then the 
following inequality must hold: 



1 



A (uj\ + Luj) + g < 0. 



Assuming for the sake of simplicity that oj\ = ic-j = ljq, we can rewrite it as 
(m 1 m 2 y 1 {m x g xx + m 2 Lolg pp ) 2 + w 2 (m 1 g px - m 2 g xp Y 



< 0. 



(4.12) 



(4-13) 



If u>q > 0, the only possibility for fulfilling this inequality is to assume that the masses have opposite 
signs. Of course, such a system is unstable, and it can be considered only as an extremely simplified 
model. Nonetheless, sometimes the models of this sort were considered [p3| . So we discuss briefly this 
case as well. Suppose for simplicity that toi = 1, ni2 = — 1. Then g < 0. The square roots in Eq. (4.7) 
can be extracted if \g\ = 2wq|A|. In this case the imaginary parts of lj± differ from zero, provided A < 0. 
Such a situation holds for the following relations between the coupling coefficients: 



9xp — 9px 5 



2 

uj o9pp- 



Then 



lj± = uj ± i-f, 
The explicit form of the matrix Ru reads 

Rn = cosh.7i 



7=|A| 1/2 



COS iOot 

ujq 1 sin wot 



/ 2 
\9px 



2 2 \V2 

^adpp) 



-ujq sin toot 

COS Uot 



(4.14) 



(4.15) 



(4.16) 



Consequently, the amplitude of the oscillations increases without bound. The drift matrix A (3.19) 
depends on time as follows: 



A = 



7 tanh jt 



1 7 tanh jt 
An interesting model described by the Hamiltonian 

H = PiP2 + LolxiX 2 + 7 (X 2 P2 - XlPl) 

was proposed for the first time by Bateman [54| . The Lagrangian 

L = ±ix 2 - {lJq + 7 2 ) xix 2 + 7 (xi± 2 - iix 2 ) 



(4.17) 



(4.18) 



(4-19) 



was considered by Morse and Feshbach |pq| . This system was investigated, e.g., in [|56| (see also flgf]). 
The equations of motion in this case read 

X! + 2 7 ii + (ljI + 7 2 ) x-i =0, x 2 + 2 7 ± 2 + (w 2 + 7 2 ) x 2 = 0. (4.20) 
Consequently, at the classical level we have damping in the first mode and amplification in the second 



mode. Moreover, at the classical level both particles are completely independent due to Eq. (4.20). The 
quantum picture is more complicated, since the quantum behaviour is governed not by the second-order 
equations of motion, but by the Hamiltonian (4.18), in which the dynamical variables of both particles 
are entangled. (The nonunique correspondence between the equations of motion and the Lagrangians or 
Hamiltonians leading to them, as well as related ambiguities of quantization, were investigated in detail 
It can be proved |57| 



in 5* 



that no Hamiltonian leading to Eq. (4.20) and coinciding with Eq. ([4.l|) 
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for 7 = exists.) In particular, the velocity of each particle is determined by the generalized momentum 
of the other particle: 

&1=P2— 7#1, X2—pi+JX 2 . 

The two other Hamilton equations read 

pi = "fPl - LUqX 2 , P2 — 1P2 - WqXi. 

In terms of the blocks of the matrix R these equations can be rewritten as follows: 



Rife aRifc + bR 2 i, 
where a and b are 2x2 matrices, 



R 



2 k 



bR 



Ik 



aR 



2k; 



7 


b = 


-ul 


-7 




1 



(fc = l,2) 



a = 

Eliminating the matrix R2fc we obtain the second-order equation 

Rife - 2aRife + ( 7 2 + wg) R u = 0. 
Seeking its solution in the form Rife. = cxp(Ai)R we obtain the characteristic equation 

A 2 - 2aA + ( 7 2 + luI) I = 0, 

whose solution reads 

A = a ± iujqL. 

Taking into account the initial conditions 

Rn(0)=I, Rn(0)=a, R 12 (0) = 0, R ia (0) = b, 
we obtain finally the matrices 



Rn = cos wot 











-7* 



R, 



sin LUgt 



-uj e 




7* 



(4.21) 



(4.22) 



Consequently, after averaging over the state of the second particle we have a system with a damped 
coordinate, but with a momentum increasing in time without bound. Moreover, the momentum will no 
longer be related to the velocity. The drift matrix (3.19) in the Fokker-Planck equation for the averaged 
Wigner function of the first particle depends on time as follows: 

ojo tan(woi) +7 

ujq tan(^oi) — 7 



We see that the reduction of "Bateman's mirror model" (4.18) does not lead to a damped quantum 
oscillator in the conventional meaning of this term. 



5 Oscillator in a thermostat: weak coupling limit 

We now proceed to a more realistic model, where the oscillator under study (its frequency will be denoted 
by ujq) is coupled to a large number of other oscillators with frequencies This model was the subject 
of investigations in many papers: see, e.g., |l^] - Pq| . A more comprehensive reference list can be found 
in the review . We assume that each oscillator is described by the standard Hamiltonian 

Hi = ^ Oi + w?ar?) 

with unit mass (this can easily be achieved by rescaling the coordinates) , while the quadratic interaction 
Hamiltonian is chosen in the most general form: 



H-,^ — 



{ziPiPo + ViPiXQ + UiXiPo + giXiXo) . 



(5.1) 
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The coupling constants z%, Vi, Ui, gi are assumed to be time independent. 

Equation (3.4) and its initial condition are equivalent to the following equations and initial conditions 
for the blocks of the matrix R of ( [3.12 ): 



(5.2) 

(5.3) 
(5.4) 
(5.5) 

The 2x2 matrix An and the 2M x 2M matrix A22 (M is the number of oscillators in the reservoir) read 



Rn 


— >4iiRii - 


- -4l2R21j 


Rn(0) 


= 1, 


R21 


= ^2lRll - 


- .A22R2I! 


R 2 i(0) 


= 0, 


R12 


= .A11R12 - 


- -4l2R22, 


Ri 2 (0) 


= 0, 


R-22 


= «42lRl2 " 


f ^422^-22) 


R 22 (0) 


= I. 



An 



-uil 




1 





-diag(o;?,...,cj?,. 



diag(l, . . . , 1, . . .) 

The matrices A12 and A21 are rectangular with the dimensions 2 x 2M and 2M x 2, respectively: 



A12 = 



Zl 



~9i 

Ui 



(5.6) 



(5.7) 



A21 = 



-U\ 


-9i 


-Ui 


-9i 


Zl 


V\ 


Zi 





(5.8) 



In this section we consider the case where all the elements of the interaction matrices A12 and A21 



are small. Then we may use perturbation theory. In the zeroth approximation we obtain from Eq. (5.2) 

R$(i) =exp(^iii). (5.9) 

Putting this expression into the right-hand side of Eq. ( [5.3] ) we obtain the first-order solution for the 
matrix R21, 



(*) = ex P (y422<) / exp (—A22T) A21 exp (Aut) dr. 
Jo 



(5.10) 



Then Eqs. (|3.19| ) and (5.2) lead to the following first-order approximation for the m atrix A governing 
the evolution of the average values of the subsystem variables according to Eq. ( 2. IS ): 



A (1) = An + A 12 R.2i (t) Rff(*) 



(5.11) 



Taking into account Eqs. (5^) and ( 5.10 ) and making the change of variable t — t = x in the integrand 
we arrive at the formula 



fj, — A — An = A12 I exp (^222;) A21 exp (— Anx) dx. 
Jo 



The solutions to Eqs. (5.4) and (p . 5|) in the same approximation read 



R 2 2 ) (t) = exp(^ 22 t), 



R12 [t) = exp (Ant) / exp (-Aut) A12 exp (^22^) dr. 
Jo 



(5.12) 

(5.13) 
(5.14) 
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Then the diffusion matrix of the Fokker-Planck equation ( 2.18 ), due to Eqs. ( 3.21 ) and (5.4), reads 
D = symfU^R^FR^t)" 



= sym |./4. 12 exp(./4. 22 t)F J exp (.A 22 t) A12 exp (-Axirj drexp^ n ; 
The explicit forms of the matrices R^(f) and R 22 (t) are as follows: 



exp (Ant) 



coswot — ujosmwot 
loq 1 sin cjo t cos ujq t 



exp (A22t) = 

If we choose the thermostat variance matrix to be 



diag(cos ujit) diag(— u>i sin ojit) 
diag(o^ sinuj it) diag(cos Uit) 



diag(a;?/i) 
diag(/0 



(5.15) 

(5.16) 
(5.17) 

(5.18) 



(in particular, F may be an equilibrium variance matrix for the thermostat variables), then it describes 
the steady-state solution of Eq. ( 2.14 ) in the absense of the interaction: 



F(i) = exp(^l 22 i)Fexp U 22 i) = F 



(5.19) 



Consequently, one may interchange the matrices exp (.4 22 t) and F in Eq. (5.15) in accordance with the 
formula 

exp (A 2 2t) F = F exp (-Avit) ■ 



ThenEq. j5~15| ) assumes the form (x = t — r) 



D = sym < / exp (Anx) A12 exp {—A22X) dx F.Ai 2 



The explicit expressions for the matrix elements of 2 x 2 matrix fi of (5.12) are as follows: 
Mi 



t ii2 = ^Y,[ WoKiS i )+w oA 4 ^ ( ) +G i c} +) 



We have introduced the notation 



g( ± ) _ sin(cj, ~UJ )t ± amjcji + Ljp)t 



LOi - UJQ 



C 



(±) 1 — COs(u>i — LO )t 1 — COs(Wj + LU )t 



UJl - UJQ OJi + UJQ 

Other parameters are the bilinear combinations of the coupling constants: 



Aj = giZi - mvi, 
Gi = ujivf + gj/u)i, 



Ki = uJiZ t Vi + giUi/uji, 
Zi = UiZi + uf/cui. 



(5.20) 
(5.21) 

(5.22) 

(5.23) 
(5.24) 
(5.25) 

(5.26) 
(5.27) 

(5.28) 
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Eq. (5.21) results in the following expressions for the elements of the diffusion matrix D: 



(5.29) 



D22 = , 



ZiS. 



(+) 



Ki C. 



(") 



(5.30) 



(-) 



(5.31) 



We see that in the general case both the drift matrix and the diffusion matrix have rather complicated 
time dependences. However, under certain conditions the formulas can be simplified, if we proceed to 
the continuum limit. This means that we assume the number of oscillators in the reservoir to be very 
large and the frequencies uji to be so close to each other that we may replace the sums by integrals over 
du>i = duj. Then we need to calculate integrals of the following type: 



sin(uj±uj )t , 
tpyuj) ■ aui, 

UJ ± UJQ 

1 — cos(w ± OJo)t . 
<p{0J) : duj. 



(5.32) 



(5.33) 



U) ± UJQ 

From physical consider ation s i t is c lear th at th e elem ents of the matrices /i and D are determined mainly 
by the terms in Eqs. ( |3.22[ )-( ^25 ) and ( 5.29|) -( 5.31 ) that correspond to frequencies near ujq, since the 
most effective interactions between the oscillator under study and the thermostat oscillators take place 
under resonance condition. Indeed, for a sufficiently smooth function ip(uj) and for t 3> uj^ 1 , only point s 
belonging to the domain \uj — ujq\ < t^ 1 make a significant contribution to the integrals ( 5.32 ) and ( 5.33| ), 
due to the rapid oscillations of the trigonometric functions outside this domain. Thus, assuming that 
lu = ujq in all the functions except sin(cj + LUo)t, we may evaluate the integral cr + as follows: 



2luq 



/ sin(cj + uJo)t doj — [cos(cji + uJo)t — cos(a->2 + <^o)t] 

Jut 2u t 



Consequently, we may neglect the value of er+ at uj^t ^ 1. The same is true for the integral 

cos{uj + uj )t (I 
<p(co) duj ~ C 

UJ + UJQ \UlQt 

As concerns the integral c_, its value does not depend on time for t ^> 0Jq , since making the substitutions 
x = u> — u>0i V = xt, we obtain 



Lp(uj ) 



sin xt 



dx — (p(too) 



sin y 

y 



dy = Tnp(uj Q ). 



(5.34) 



To evaluate the integral J (/3(a;) cos ^~"°'* t du> we use the Taylor expansion cp(u) = f(uJo) + (p'(ujq)(uj — 
luq) + ■ ■ ■, make the substitution x — to — ujq, and expand the limits of integration from —00 to 00. Then 
the first integral vanishes, because the function cos(xt)/x is odd. The second integral decreases at least 
as l/t when t — > 00. Therefore, we may assume that the integrals S± do not depend on time at Loot 1: 



UJ — UJQ 



duj. 



(5.35) 



In the case of the "minus" sign the principal value of the integral is implied: it is designated with the 
symbol J . We arrive at the following expressions for the matrix elements of the matrices \x and D in the 
continuum limit: 



Mil 



A* 12 



1 / n r a / \ 1 ^, Nl f' UK,(uj)v(uj) . 

--m{ojo)[A{uj )+ujv 1 G(ujo)]+ — : \ duj, 

Z I UJ — UJn 



-■nv(uj )uj Q n{uj ) 



v(uj) 



[woA(w) +ujG(uj)] duj, 



(5.36) 
(5.37) 
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1 f f v(oj} 
M21 = ^v{u)o)lJq X k{ljq) - / —5 2 [AM + w ^( w )] dw, (5.38) 

l J LU — UJq 

M22 = -^1/(0/0) [AM + ^M] - /' dw, (5.39) 



2 



On = J7ri/(wo)/(wo)wo [woA(wo) + G(wo)] - ^0 / ^M^^ (5 4Q) 



2 



"0 



^22 = iWoJ/M [A(w ) + uj Z(uj )} + f ^P^M duJi (5 . 41) 
D 12 = - ^(woM/Mk^o) + [ wf ^ v ^ [u? Q Z{u) - G(w)] dcj. (5.42) 



2 

Here f(tL>) is the density of states function, whereas the functions Z(u>), G(cj), etc. are obvious general 



izations of functions defined in Eq. (5.28) 



In principle, the frequency dependences of the coupling constants can be chosen in such a way that 
the integrals in Eqs. ( 5.36| )-( |5.42j ) vanish. For instance, this is possible provided the corresponding 



combinations of functions K, v, /, G, Z , A/lu, understood as functions of the argument x — u , do not 
change their values under the reflection in the point xq — ujq, and these functions (or at least the density 
of states) decrease sufficiently rapidly with distance from the point xo in both directions. In such a case 
all the coefficients p,ik and Dik are determined by the values of the aforementioned functions at the point 
ujq. Moreover, the diffusion coefficients are proportional to the corresponding elements of the matrix p: 

Du = -o>o/oMhj D 22 = -foV22, Dx2 = -/0M12 = -Wq/oM21' (5-43) 
We see that in the continuous weak coupling limit the reduced Wigner function of the oscillator obeys 



(under certain conditions) the Fokker-Planck equation (2.18) with time-independent coefficients at times 
t 3> w^ 1 - Let us check, however, whether condition ( [2.27 ) is fulfilled. Since TrA = Tr/i, then due to Eq. 



(5.43) we must check the inequality 

w o/oVnM22 - /0M12 ^ ^ (Mn + A^2) 2 /16. (5.44) 



Taking into account Eqs. ( 5.36| )-(5.39) and Eq. (5.28), we arrive at the inequality (all the functions are 
taken at the point u) = ljq) 

(Afuj/h) 2 [2A 2 + A(wZ + G/lu)] > [2A + (uZ + G/cu)} 2 . (5.45) 

It cannot be satisfied for an arbitrary choice of coupling constants. For instance, it is violated if 
A = 0, or when any three of four coefficients g, 2, u, v vanish. 

This result seems paradoxical. Indeed, we started from the exact equation of motion for the density 
matrix of a closed system, found the exact solution to this equation, and after this we performed averaging 
over the thermostat degrees of freedom. Since the laws of quantum mechanics were not violated at any 
step, the reduced density matrix must be positive definite at any time for quite arbit rary coupling 
constants. On the other hand, if, for instance, Zj, = Ui = Vj = 0, then inequality ( 5.45| ) does not hold. 



Hence, following the reasonings given in Sec. |2j, we could obtain a nonpositivc definite density matrix in 
the process of evolution! 



This apparent contradiction is resolved in the following way. Inequality (2.27) is a necessary and 
sufficient condition ensuring that any density matrix that was positive definite at any instant of time 
will remain positive definite at all subsequent moments. But in the problem under study we have the 
selected instant t — 0: this is just the moment when the interaction with the thermostat was turned on. 
Since in the presence of the interaction with the environment the evolution of the oscillator density matrix 
is nonunitary, the set of density matrices p(t) arising from all initially admissible density matrices p(0) 
does not coincide with the set of all admissible density matrices. In particular, correct initial density 
matrices cannot turn into the specific ones that could become nonpositive definite at some instant of 
time in the case of violation of inequality ( 5. 45] ). Therefore, there is no need to check conditions like 



(2.27) or (5.44), (5.45) when the density matrix of the subsystem is obtained by reduction of the exact 
density matrix of the closed system: the reduced density matrix (calculated with the proper accuracy) 
turns out to be positive definite automatically. 
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However, if the goal is the derivation of a self-consistent Fokker-Planck equation on the basis of a 
"microscopic" model of the oscillator interacting with a large reservoir, then we have to recognize that 
the underlying "microscopic" model cannot be quite arbitrary: its parameters must satisfy rather strong 
restriction (5.45), in order to prevent the appearance of the unphysical solutions when this equation is 
applied to arbitrary initial states. 

Returning to the analysis of Eqs. ( 5.36| )-(5.43), we notice that the matrix 



wg/o 
fo 



satisfies, due to Eq. ( |5.43 ), the relation 



AF + F A + 2D = 0. 



(5.46) 



(5.47) 



This means that Fo is the steady state solution to Eq. (2.19) for the oscillator variance matrix, inde- 
pendently of the concrete values of the coefficients of the drift matrix fi^. In particular, if 



/o = ^cothfe 



2ui 



\2kT 



(5.48) 



then the matrix (5.46) coincides with the equilibrium variance matrix of the oscillator. Consequently, the 
steady state solution of the Fokker-Planck equation with the coefficients ( 5.36| )-(5.42) is the equilibrium 
Wigner function of the oscillator, if the thermostat itself is in the equilibrium state. Moreover, we may 
imagine a situation with 



fo = 7, — ai coth 
2u ^ 



( 



flLOn 



\2kT t 



a,, > 0, 



(5.49) 



Then, the oscillator under study exhibits relaxation to the Gaussian steady state with the variance matrix 
given by Eqs. (5.46) and (5.49). Such a situation may be realized when the thermostat consists of several 
large independent subsystems possessing their own temperatures (see, e.g., [[59|). In this case n is the 
number of subsystems, and on is the "weight" of each subsystem. 

We see that the final steady state of the oscillator weakly interacting with a thermostat does not 
depend on the concrete values of the coupling constants, provided the thermostat was initially in an 
equilibrium or quasi-equilibrium (described by Eq. ( 5.49| ) state. Now let us analyze possible forms of the 
drift matrix A. Due to the property — W5V21 its characteristic equation det(A — AI) = reads 



A 2 - (/in + M22) A + uj 2 + /in/122 - M12M21 = 
The solutions to this equations can be written as Ai,2 = —7 ± i^*, with 

1 



0. 



1 

"2 TrM 



w o - 7 _ M22) - M12M21 



\ (Mil + M22) = ^7w(w ) [ 2z ^o + wqZo + ^0 lG o] 



(5.50) 



(5.51) 



(5.52) 



The last formula holds even when the integrals in Eqs. ( 5.36| ) and (5.39) are not equal to zero. Assuming 
the damping coefficient 7 to be positive and comparing Eqs. ( 5.52| ) and (5.45), we may rewrite inequality 
(5.45) as follows (recall that all functions are taken at u> — luq): 

{4fuj/h) 2 A > 2A + (ujZ + G/uo) = 4 7 /[7n/(w)]. (5.53) 

In particular, for an equilibrium thermostat at zero temperature (when / = H/2lu) we obtain the relation 

2A> luZ + G/w, (5.54) 



which due to Eq. ( 5.28 ) turns into the inequality 

(u + v f + (ujz - g/uo) 2 < 0. 



(5.55) 
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Thus we arrive at a striking conclusion: the microscopic model of the interaction between the oscillator 
and the equilibrium reservoir at zero temperature, based on the general quadratic Hamiltonian (5.1), leads 
to the Fokker-Planck equation with time independent coefficients, valid for any physically admissible 
initial states of the oscillator, if and only if the coupling constants at the main oscillator frequency satisfy 
the restriction 



u 



-v , z = go/col 



In this case Eqs. ( 5.36 )-( 5~39|) yield 



Mil = M22, M12 = M21 = 0, 



so we have the unique drift matrix 



(5.56) 



(5.57) 



(5.58) 



with the damping coefficient 



7 = 7w(w )Ao, A = ul + g\ju\. 



(5.59) 



Due to Eq. (5.57) the frequency oj* in formula (5.51) exactly equals the oscillator eigenfrequency ujq. 
The equations of motion for the average values of the coordinate and momentum read 



-IP-UqX, (5.60) 

x = p — jx. (5.61) 

One can see that consistent quantum mechanical consideration do not result in the conventional classical 
equations ( 2.16Q . To understand the origin of Eqs. ( |5.60 ) and ( |5.61 ), let us introduce the annihilation 
and creation operators 



UJqX + Ip 



>+ UJqX — ip 
1 I — 



(5.62) 



y/2hu)Q \f2TlLOQ 

which are the most natural for the description of a quantum oscillator. It turns out that precisely Eqs. 



(5.60), (5-61) lead to uncoupled equations for a and a' 



-iloqcl — 7a, 



ILUqCL' 



7a 1 



(5.63) 



Denoting the annihilation and creation operators for the thermostat oscillators by Cj and cj we may 
rewrite the interaction Hamiltonian (|5.l[ ) (for thermostat oscillators whose frequencies coincide with ojo) 
as follows: 



E 



SqSl'c, 



(5.64) 



5 a = ga/uj Q + iuQ, A = \S \^ 



Hamiltonian (5.64) was considered in almost every paper devoted to the models of a quantum damped 
oscillator (see, e.g., review Q). But frequently it was chosen only because it is the simplest one. We 
have shown in fact that it is the only possible quadratic interaction Hamiltonian ensuring the validity of 
the resulting Fokker-Planck equation for any initial states taken at any initial times. This is probably 
related to the quantum nature of the interactions between the systems: each act of interaction must 
consist in the annihilation of a quantum in one system and its creation in another system. Precisely the 



Hamiltonian (5.64) expresses this property in the most distinct form. 



Nonetheless we have no sufficient grounds for claiming that other interaction (quadratic) Hamiltonians 
should be excluded. But they will result in the Fokker-Planck equation with time-dependent drift and 
diffusion matrices describing nonexponential relaxation. 

Recall that we assumed the integrals in Eqs. ( 5.36 )-( 5.42 ) to equ al ze ro. What will happen if 
we abandon this assumption, but assume instead that the constraints ( |5.56j ) hold for all the coupling 
constants? In this case we have k(u>) = at all frequencies. Furthermore, G(u>) = u> 2 Z(lj) = wA(w). 
Consequently, the integral terms "survive" only in the off-diagonal elements of both matrices fj, and D: 



A* 12 



Mil = M22 = 



du>, 
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<■' KoQAM A 

M21 = - / —5 — dw, 



W0-D22, 



Do 



7/(^o), 





D 12 



(5.65) 



(the function 7(0;) is defined by Eq. ( 5.59 ) with loq replaced by oj). 

The presence of nonzero coefficients /112 and /i2i may be interpreted as some kind of renormalization 
of mass and eigcnfrcqucncy of the main oscillator due to the interaction with the environment. However, 
inequality (2.26) is obviously violated at zero temperature, when / = H/2lu, since due to Eq. (5.65) 
coefficient D12 is strictly negative for all temperatures (evidently, both functions v(lo) and f(uj) are 
positive). This example shows once more that the self-consistent Fokker-Planck equations with time- 
independent coefficients can be derived from microscopic models only in exceptional cases. 



6 Oscillator in a magnetic field. Weak coupling with a thermo- 
stat 

We now consider, within the framework of the same scheme, the case where the subsystem under study is 
a two-dimensional isotropic oscillator with eigenfrequency uiq and mass m placed in a uniform magnetic 
field TL characterized by the cyclotron frequency 



The Hamiltonian of this subsystem reads 

1 



Ho = 



2m 



= eH/i 



(tt 2 x + % 2 y ) + -mwg (x 2 + y 2 ) 



(6.1) 



(6.2) 



where n x and tt v are the operators of the kinetic momentum projections, related to the canonical mo- 
mentum p and the vector potential A in the usual way: 



(e/c)A(x,y), 



ihmuj. 



(6.3) 



The problem of constructing the Fokker-Planck equation for this subsystem is reduced to that solved 
in the previous section, because the Hamiltonian ([T^) can be expressed as a sum of two oscillator 
Hamiltonians: 



1 



1 



The annihilation operators 



H = hcu+ [a 1 a+ - + [b [ b+ - 



a = (2rnfin) 1 ^ 2 [tt x + ijr y + moj- (y — ix)] , 
b = (2mhil)^ 1 ^ 2 [tt x — iiTy — muj+ (y + ix)] 



satisfy the commutation relations 

The frequencies are defined as follows: 

u)± = — (fi ± oj), ft 



a, b 



2\l/2 



(6.4) 

(6.5) 
(6.6) 

(6.7) 
(6.8) 



In the continuous weak coupling limit the self-cons isten t equations of motion for the first-order average 
values of the operators a and b are given by Eq. ( 5.63 ), provided one replaces loq with lo + and 
Furthermore, two different damping coefficients are possi ble: t hey a re det ermined by the density of states 
and coupling constants at the frequencies uj± (see Eqs. (|5.59|) and ( 5.64| )): 



7± = kv{!jJ±)\5(uj±){' 



(6.9) 
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The relations inverse to Eqs. (6.5) and ([T(]) read 

it x = {mh/2tt) 1/2 uj + (a + a f ) + uj- (b + S f J 

uj+ (a) - a) + lu- (b - 



f y = i{mh/2VL) 1/2 

x = i(h/2mn) 1/2 
y = (h/2mfl) 1/2 



a - a f + b - P 
a + a) -b-tf 



(6.10) 

(6.11) 
(6.12) 
(6.13) 



The average values of the coordinates and the kinetic momenta obey equations resulting from equations 
of the form of ( 5.63| ): 

tt x — —an x + uiTTy — mui 2 x + mu>Qey, (6.14) 

(6.15) 
(6.16) 
(6.17) 



TT y = —ujit x — aiTy — muj 2 ex — mujQy, 



x ■ 

y 



m w x 



m en v 



m eir x 



T)X, 

ryy. 



(7_- 7+ )/n. 



We have introduced the notation 

a = (7+w + + 7-cj_) /O, ?/ = (7+w_ + 7_w + ) /O, 
The second-order equations of motion read 

x + (7- + -ujy+ (u/q + 7-7+) # - (7-^+ - 7+W-) y = 0, 
2/ + (7- + 7+) 2/ + + (cJq + 7-7+) y + (7- w + - 7+ w -) x = 0. 



(6.18) 

(6.19) 
(6.20) 



We see that "one-photon" interaction with a thermostat of the form of (5.64) results in coordinate- 
dependent forces perpendicular to the vector r = (x, y) and proportional to the damping coefficients. 
The necessity of introducing such forces was shown earlier in [[39], ^| within the framework of a 
phenomenological approach. As was demonstrated in these papers, if the second-order equations of 
motion (for a charged particle or oscillator placed in a uniform magnetic field) contain the term —71", 
then it is impossible to satisfy simultaneously Eq. (2. IE) with the equilibrium matrix M(T) for all 
temperatures (including T = 0) and condition (2.20), unless a force of the form f = [h x r] is introduced. 
Now we have arrived at the same result on the base of the microscopic model. Moreover, the relation 
between the velocity and the kinetic momentum becomes much more complicated than in the conservative 
case (see Eqs. (6.16) and (6.17). 

It is clear from the preceding section that the steady state Wigner function at t — ► 00 coincides with 
the equilibrium distribution that was found in refs. j^, [Tq |. Since this distribution is Gaussian, it is 
completely determined (see Eq. (3.8)) by the equilibrium variance matrix 



M« 



Mir 








M a 







M* 


-M a 










-M 


a M p 







M a 








Mp 












(- 


-) 


sinh f2 — 2 


— sinh 


,7; 




n 2 




VL 



M p = 



h sinh 57 
mflQ ' 



Ma= 2Q 



— sinh fi — sinh Co 



where 



Q(/3) = cosh £1 — cosh a), 
and (3 = 1/kT is the inverse temperature of the thermostat 



-Phu>, 



(6.21) 

(6.22) 

(6.23) 
(6.24) 

(6.25) 
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The drift matrix A corresponding to Eqs. 



([[l|)-([[rf) reads 



-Q 
-LO 



LO 

—a 



-TTVjJ\ 



-V 



(6.26) 



Putting this matrix into Eq. (5.47) with the matrix A4^ e ^ instead of Fp we obtain the diffusion matrix 



D 
















D v 


-D a 








-D a 















Dp 



(6.27) 



with the following coefficients, 

mh 



Ar = 



2^g 

2fiQ 



p 2mnQ 



(7+W+ + 7_cj 2 _) sinhfi — (7+cJ^. — 7_lj 2 _) sinhu) 
(7+W+ — 7_oj_) sinhfJ — (7 + cj + + 7_lj_) sinhcj 

h 



(7+ +7-) sinhfi — (7+ — 7_) sinhcD 



(6.28) 

(6.29) 
(6.30) 



In particular, at zero temperature (/3 = 00) we obtain 

£>(low) 



2^(7+^ +7-W-), 



M low > = 4(7+"+ -7-"-), 



(low) _ 



2ft 



(7+ +7-) ■ 



2mfi 

In the opposite, high-temperature, case (/? — > 0) we have 



2 1 



(high) 



-kTe. 



(6.31) 



(6.32) 



It is noteworthy that all three diffusion coefficients remain nonzero even in the high temperature limit, 
which is usually identified with the quasiclassical limit. Recall that in classical statistical mechanics it is 
usually implied that the only nonzero diffusion coefficient is D^. 

Various sets of the diffusion coefficients compatible with inequality ( 2.20| ) and leading to an equilibrium 
steady state (with the variance matrix (6.21)) in the limit of infinitely small damping were constructed 
within the framework of the phcnomcnological approach in (3^]. However, none of them had a structure 
similar to that given by Eqs. j6^)-(|63c|). 

For example, only the coefficient D^ was proportional to the temperature in the high-temperature 
limit, whereas the other diffusion coefficients decreased as 1/kT, in contrast to Eq. (3.32). This difference 
is due to at least two causes. First, it was assumed in 39[ 48[ |i~9| that the elements A41 and A32 of the drift 



matrix must be zero, although the elements An and A23 could be nonzero. Eq. (6.26) shows that within 
the framework of the microscopic approach such a choice is impossible, since all these coefficients are 
proportional to the parameter e. Furthermore, in the aforementioned papers we admitted the possibility 
that some coefficients of the drift matrix (those related to the damping) could depend on temperature. 
In principle, such a possibility (i.e., the time dependence of the coupling constants in the interaction 
Hamiltonian) is not excluded within the framework of the microscopic approach as well. Then the high 
temperature limit of the diffusion matrix coefficients could be quite different from Eq. (3.32). 

The expressions (6.26)-( 6.30| ) are simplified in the special case of the equal damping coefficients, when 
7+ = 7- = 7o: 



A (o) 



-70 


LO 


—mLOQ 


— LO 


-70 










-70 





m _1 








— IJILO, 



-70 



(6.33) 
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!><?> = 



2Q.Q 



2lOq) sinhfi 



cjQ sinh a) 



^ 0) = 



^7o 
20Q 



) sinh fi — f2 sinh c 



2,(0) 



^70 



sinh £7. 



m£2Q 

If coefficient 7_ tends to zero sufficiently rapidly as w_ 



(6.34) 
(6.35) 
(6.36) 



0, then the set of diffusion coefficients 
(6.28)-(6.30) possesses the finite limit for a free particle in a magnetic field, when luq = ld_ = 0, fi = u>, 
7+(w+ = w)=7: 

D„ = — ^mTiij cotho), (6.37) 



D„ 



D„ 



1 



7?icoth(D, 



7?i 



a 



7) 



2mcj 
»7 = 0, 



cotha), 



-7/0;. 



(6.38) 

(6.39) 
(6.40) 



In this case the operators b and w become the integrals of motion, whose real and imaginary parts 
are connected with the center- of- orbit operators in a uniform magnetic field [p0| -[|62|. 

One should remember, however, that the results obtained in this section can be justified only under 
rather strong limitations imposed on the interaction Hamiltonian. First, it must be written in the specific 
form (5.64) at the resonant frequencies. Secondly, the off-resonance terms must ensure the disappearance 
of the integral terms in Eqs. (5.36)-(5.42). In particular, the density of states must decrease sufficiently 
rapidly as lu — > oo. 

All these conditions are violated, for example, in the case where the role of a reservoir is played by a 
quantized electromagnetic field coupled to the oscillator by means of the standard interaction Hamiltonian 
in the dipole approximation: 



H: 



int = -^(2^) 1/2 E% 1 [a(k^)+ct(k„a)] 



(6.41) 



where c(kj, a) is the operator of annihilation of a photon with wave vector kj, frequency de- 
polarization fi; Tj.ct is the unit polarization vector perpendicular to the vector k, . The density of states 



c|kj|, and 



is proportional to u> 2 in this case, and the integrals in Eqs. (5.36)-(5.42) diverge. Consequently, in this 
case the radiation damping leads to nonexponential relaxation. 
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